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ABSTRACT 

We present a study of the dynamical properties of 125 compact stellar systems (CSSs) in the nearby giant 
elliptical galaxy NGC 5128, using high-resolution spectra (R « 26 000) obtained with VLT/FLAMES. Our 
results provide evidence for a new type of star cluster, based on the CSS dynamical mass scaling relations. All 
radial velocity ( v r ) and line-of-sight velocity dispersion (cti os ) measurements are performed with the penalized 
pixel fitting ( ppxf) technique, which provided <r ppx f estimates for 115 targets. The <r ppx f estimates are corrected 
to the 2D projected half-light radii, <ti/ 2 , as well as the cluster cores, <To, accounting for observational/aperture 
effects and are combined with structural parameters, from high spatial resolution imaging, in order to derive 
total dynamical masses (Addyn) for 112 members of NGC5128’s star cluster system. In total, 89 CSSs have 
dynamical masses measured for the first time along with the corresponding dynamical mass-to-light ratios 
(T p "). We find two distinct sequences in the T p ) n -Add yn plane, which are well approximated by power 
laws of the forms Yy yn oc Ad° y ^ ±004 and Yy yn <x Ad° y 7 ® ±0 ' 04 . The shallower sequence corresponds to the 
very bright tail of the globular cluster luminosity function (GCLF), while the steeper relation appears to be 
populated by a distinct group of objects which require significant dark gravitating components such as central 
massive black holes and/or exotically concentrated dark matter distributions. This result would suggest that 
the formation and evolution of these CSSs are markedly different from the “classical” globular clusters in 
NGC 5128 and the Local Group, despite the fact that these clusters have luminosities similar to the GCLF turn¬ 
over magnitude. We include a thorough discussion of myriad factors potentially influencing our measurements. 
Subject headings: galaxies: individual: (NGC 5128) - galaxies: star clusters: general - galaxies: spectroscopy 
- galaxies: photometry 


1. INTRODUCTION 

Globular clusters (GCs) are among the oldest stellar sys¬ 
tems in the Universe (Krauss & Chaboyer 2003). They have 
witnessed the earliest stages of star formation and were also 
present during later epochs of structure formation. Apart from 
resolved stellar population studies of galaxies, which are re¬ 
stricted primarily to the Local Group, extragalactic globular 
cluster systems (GCSs) provide one of the best probes to 
investigate the formation and assembly histories of galaxies 
(Harris 1991; Ashman & Zepf 1998, 2008; Peng et al. 2008; 
Georgiev et al. 2010). Various avenues of study can be em¬ 
ployed to this effect, including the analysis of GCS kinemat¬ 
ics, their metallicity distribution functions, chemical enrich¬ 
ment histories, age spreads, and/or combinations thereof. 

At a distance of 3.8±0.1Mpc, (Harris et al. 2010), corre¬ 
sponding to an angular scale of 18.5pc arcsec -1 , NGC5128 
(a.k.a. Centarus A) is the nearest giant elliptical (gE) galaxy 
to the Milky Way (MW), yet it is still too far for large scale 
star-by-star investigations to be technologically feasible. For¬ 
tunately, much has been learned about this galaxy from its 
rich GC system. 564 of its GCs have radial velocity confir¬ 
mations, and others have been confirmed, for example, via 
resolution into individual stars (van den Bergh et al. 1981; 
Hesser et al. 1984, 1986; Harris et al. 1992; Peng et al. 2004; 
Woodley et al. 2005; Rejkuba et al. 2007; Beasley et al. 2008; 

t Based on observations collected under program 081 .D-0651 (PI: Ma- 
tias Gomez) with FLAMES at the Very Large Telescope of the Paranal Ob¬ 
servatory in Chile, operated by the European Southern Observatory (ESO). 
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Woodley et al. 2010a). This GC sample, thus, rivals the en¬ 
tire population of GCs harbored by the Local Group, despite 
there being ~ 600 —1400 GCs still left to find/confirm in the 
halo regions of NGC 5128 (Harris et al. 1984, 2002b, 2010, 
2012). Notwithstanding this incompleteness, previous studies 
have already shed much light on the photometric, chemical, 
and kinematical properties, as well as the past and recent for¬ 
mation history of this massive nearby neighbour. 

GCs are well known to inhabit a narrow range of space de¬ 
fined by structural parameters such as half-light, tidal and core 
radii ( rh , r t and r c , respectively), concentration parameter 
c = log ( rt/r c ) (King 1966), velocity dispersion (a), mass-to- 
light ratio (T), etc. called the “fundamental plane” (Djorgov- 
ski 1995). Studies of Local Group GCs (e.g. Fusi Pecci et al. 
1994; Djorgovski et al. 1997; Holland et al. 1997; McLaugh¬ 
lin 2000; Barmby et al. 2002, 2007) have shown that at the 
high-mass end of the fundamental plane, peculiar GCs such as 
wCen and Gl, the largest GCs in the MW and M31, begin to 
emerge. For example, both of these GCs show significant star- 
to-star [Fe/H] variations and are among the most flattened of 
Local Group GCs (White & Shawl 1987; Norris & Da Costa 
1995; Meylan et al. 2001; Pancino et al. 2002) and at least in 
the case of wCen, harbor multiple stellar populations with an 
extended chemical enrichment history (Piotto 2008a,b). 

Unlike the Local Group, the sheer size of NGC5128’s GCS 
generously samples the high-mass tail (> 10 6 Mg) of the 
globular cluster mass function (GCMF, Harris et al. 1984, 
2002a; Martini & Ho 2004; Rejkuba et al. 2007; Taylor et al. 
2010). Due to their intense luminosities, these massive GCs 
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are very accessible observationally, and thus provide excellent 
probes to study the formation history of NGC5128. Many of 
the most massive NGC5128 GCs show a more rapid chemi¬ 
cal enrichment history than Local Group GCs (Colucci et al. 
2013), and exhibit significantly elevated dynamical mass-to- 
light ratios (Ty yn ) above dynamical masses, A4d yn — 2 • 10 6 
M 0 (Taylor et al. 2010). This sharp upturn of Yy yn is consis¬ 
tent with a trend found by Hasegan et al. (2005) and Mieske 
et al. (2006, 2008a) in other extragalactic GCSs, and re¬ 
quires either non-equilibrium dynamical states, such as rota¬ 
tion or pre-relaxation (Varri & Berlin 2012; Bianchini et al. 

2013) , younger than expected stellar components (e.g. Be- 
din et al. 2004; Piotto 2008a), exotic top- or bottom-heavy 
stellar initial mass functions (IMFs; e.g. Dabringhausen et al. 
2008, 2009; Mieske et al. 2008b) or/and a significant contri¬ 
bution by non-baryonic matter or massive central black holes 
(BHs). While there is an ongoing debate whether the latter 
two options are valid for Milky Way GCs (see e.g. Conroy et 
al. 2011; Ibata et al. 2013; Liitzgendorf et al. 2011; Strader 
et al. 2012; Lanzoni et al. 2013; Sun et al. 2013; Kruijssen 
& Liitzgendorf 2013), there is growing evidence for the pres¬ 
ence of 10 5 — 10 8 Mr., BHs significantly affecting the dynam¬ 
ics of similarly structured, albeit more massive, ultra-compact 
dwarf galaxies (UCDs; Mieske et al. 2013; Seth et al. 2014). 

NGC5128’s GCS has been shown to follow trends sim¬ 
ilar to other giant galaxies. In particular, it has a multi¬ 
modal distribution in color and metallicity (e.g. Harris et al. 
2002b; Peng et al. 2004; Beasley et al. 2008; Woodley et al. 
2010b), corresponding to at least two and possibly three dis¬ 
tinct GC generations. Moreover, the prominent dust-lane and 
faint shells in the galaxy surface brightness distribution (Ma- 
lin et al. 1983), along with a young tidal stream (Peng et al. 
2002) provide significant evidence for recent merger activity 
on a kpc scale, while on smaller scales indications of strong 
tidal forces are seen in the form of extra-tidal light associated 
with individual GCs (Harris et al. 2002a). 

Recent models support the notion that the bulk of the star 
formation leading to massive elliptical galaxies is complete by 
z~ 3 (i.e. the first few Gyr of cosmic history), while it takes 
until z < 0.4 before ~ 80% of the mass is locked up after the 
accretion of as many as five massive progenitors (e.g. De Lu¬ 
cia et al. 2006; De Lucia & Blaizot 2007; Marchesini et al. 

2014) . Studies based on Hubble Space Telescope (HST) data 
in the mid- to outer-halo regions of NGC5128 generally con¬ 
cur with this view, in that the majority (~ 80%) of the stellar 
population is ancient (> 11—12 Gyr) and formed very rapidly, 
as evidenced by [a/Fe] ratios approaching or exceeding twice 
solar values (e.g. Harris et al. 1999; Harris & Harris 2000, 
2002; Rejkuba et al. 2011). This older population is comple¬ 
mented by a significantly younger component, forming on the 
order of a few Gyr ago (e.g. Soria et al. 1996; Marleau et al. 
2000; Rejkuba et al. 2003). 

In this paper we use velocity dispersion estimates based on 
high-resolution spectra to derive dynamical masses for a large 
sample of NGC 5128’s GCS (see e.g., Chilingarian et al. 2011, 
who carried out similar work on ultra-compact dwarfs in the 
Fornax cluster). We combine the newly derived dynamical 
information with well-known luminosities from the literature 
to probe the baryonic makeup and dynamical configurations 
of the CSSs. The results are then used to classify several dis¬ 
tinct CSS/GC populations, which are discussed in the context 
of likely origins, with potential consequences for GCSs that 
surround other gE galaxies. 


This paper is organized as follows. § 2 describes the obser¬ 
vations made as well as an outline of the data reduction steps 
taken to produce high-quality spectra. § 3 contains informa¬ 
tion on the analysis that was undertaken on the new spectro¬ 
scopic observations, as well as structural parameter data from 
the literature with which our new measurements were com¬ 
bined. § 4 discusses our results by using sizes, masses and 
mass-to-light ratios of GCs/CSSs to develop several hypothe¬ 
ses on the origins of the various cluster sub-populations that 
we find. The main text concludes with § 5, which summarizes 
our new measurements and results. Following the main text, 
we present in the appendix multiple detailed tests which rule 
out spurious results due to several possible sources includ¬ 
ing poor data quality, data analysis biases, fore/background 
contamination, target confusion, and others. We adopt the 
NGC5128 distance modulus of (m. — M) o = 27.88 ± 0.05 
mag, corresponding to a distance of 3.8 ± 0.1 Mpc (Harris 
et al. 2010), as well as the homogenized GC identification 
scheme of Woodley et al. (2007) throughout this work. 

2. OBSERVATIONS 

During five nights in lune/July 2008, 123 of the brightest 
GCs around NGC 5128 were observed using the Fibre Large 
Array Multi-Element Spectrograph (FLAMES) instrument at 
the Very Large Telescope (VLT) on Cerro Paranal, Chile. 
FLAMES is a multi-object spectrograph mounted at the Nas¬ 
myth A focus of UT2 (Kueyen). The instrument features 132 
fibres, each with apertures of 1.2" diameter, linked to the 
intermediate-high resolution (7000 <R< 30 000) GIRAFFE 
spectrograph, with an additional eight 1.0" fibres connected to 
the high resolution (R ss 47 000) UVES spectrograph mounted 
at the Nasmyth B focus, thus, allowing for simultaneous ob¬ 
servations of 139 targets 3 over a 25' diameter field of view. 

2.1. Instrumental Setup 

FLAMES often suffers from one or more broken fibres, and 
these observations were no exception. The first night of ob¬ 
servations, comprising the first three of eleven 2 400 second 
long observing blocks (OBs), were conducted with the first 
of two FLAMES fibre positioner plates, which suffered from 
two broken fibres, while the last eight OBs used the second 
plate, with only a single, but distinct, broken fibre. For this 
reason, the total number of targets observed was 138, includ¬ 
ing eight fibres fed to UVES, and 130 to GIRAFFE. Unfor¬ 
tunately the UVES targets were of too low data quality for 
useful measurements to be derived, and so we do not include 
them in the present work. From the GIRAFFE fibre budget, 
13 were allocated to recording the sky contribution to the tar¬ 
get signals. These sky fibres were used near the end of the data 
reduction process to perform an adaptive sky subtraction. In 
summary, we obtained spectra of 117 bright GCs using GI¬ 
RAFFE in the high resolution (II ~ 25 900) mode, recording a 
single AA = 21.3 nm echelle order centred at A = 525.8 nm, 
covering the wavelength range 515.2 < A < 536.5 nm. 

Table 1 summarizes our observations. The overall exposure 
times were 26 400 s; however, three GCs (GC 0310, GC 0316 
and GC 426) inevitably suffered from broken fibres, limiting 
the total integration times to 19 200 s, 7 200 s, and 7 200 s, re¬ 
spectively. In the case of GC 0058 a single exposure had to 
be omitted due to the unfortunate coincidence of a significant 

1 In principle 132 GIRAFFE fibres are allocatable, but only 131 are fully 
covered on the detector; therefore, the sum of GIRAFFE+UVES available 
fibres is 139. 
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Fig. 1.— Locations of the GIRAFFE fibre placements in relation to the host galaxy, NGC 5128. Grey points show the placement of fibres for “classical” GCs, 
while blue pentagons and red squares illustrate the positions of DGTO and DSC GCs, respectively (see Figure 6 and corresponding text in Sect. 4.1 for details 
on their definition). Green points mark the positions of sky fibres. The light red shaded region shows the approximate extent of the 25' FLAMES field-of-view, 
corresponding to ~ 27.5 kpc at the distance of NGC 5128. 


detector defect lying directly in the middle of the Mg/; triplet, 
hence the total integration time for that GC was limited to 
24000 s. 

The on-sky locations of the GIRAFFE fibres are indica¬ 
ted in Figure 1, over-plotted on an archival DSS image 4 of 
NGC 5128. Overall, the observing conditions were good for 
this observing program: The images were taken at airmass 
values ranging between 1.054 and 1.635, under seeing condi¬ 
tions in the range 0.48" to 1.42". For the purpose of correcting 
our line-of-sight velocity dispersion measurements, based on 
our final stacked spectra, for aperture and observational ef¬ 
fects (see § 3.5), we use the mean seeing value from all 11 

4 Based on photographic data obtained using The UK Schmidt Telescope. 
The UK Schmidt Telescope was operated by the Royal Observatory Edin¬ 
burgh, with funding from the UK Science and Engineering Research Council, 
until 19S8 June, and thereafter by the Anglo-Australian Observatory. Origi¬ 
nal plate material is copyright © the Royal Observatory Edinburgh and the 
Anglo-Australian Observatory. The plates were processed into the present 
compressed digital form with their permission. The Digitized Sky Survey 
was produced at the Space Telescope Science Institute under US Government 
grant NAG W-2166. 


OBs of 0.85", since all targets were observed simultaneously 
under identical seeing conditions. 

2.2. Basic Data Reduction and Cleaning 

The basic data reduction steps (bias subtraction, flat field¬ 
ing, and wavelength calibration) were carried out by the GI¬ 
RAFFE pipeline 5 . Separate calibration frame sets were used 
for each of the five nights. The pipeline recipe masterbias 
created the master bias frame from an average of five indi¬ 
vidual frames and masterflcit produced the master flat from 
an average of three bias-subtracted Hatfields. The fibre local¬ 
izations were visually confirmed to be accurate to within 0.5 
pixels for each of the 11 OBs, less than the suggested 1 pixel 
maximum to ensure accuracy. The recipe giwavecalibration 
derived the wavelength calibrations. For all five calibration 
sets, it was necessary to edit the slit geometry tables in or¬ 
der to eliminate ‘jumps’ in the final, re-binned, wavelength 
calibrated arc-lamp spectra - an extra step that is not uncom- 

5 http://www.eso.org/sci/software/pipelines 
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TABLE 1 

Star Cluster Observations 


ID 

a 

[J2000] 

5 

[J2000] 

R 

[mag] 

Vo 

[mag] 

Exp. Time 
[s] 

S/N 

GC 0028 

13 24 28.429 

-42 57 52.96 

19.65 

19.80±0.01 

26400 

2.87 

GC 0031 

13 24 29.700 

-43 02 06.43 

19.55 

19.75±0.01 

26400 

1.88 

GC 0048 

13 24 43.586 

-42 53 07.22 

19.33 

19.51T0.01 

26400 

1.79 

GC 0050 

13 24 44.575 

-43 02 47.26 

18.90 

18.74±0.01 

26400 

3.75 

GC 0052 

13 24 45.330 

-42 59 33.47 

18.91 

18.98±0.01 

26400 

5.55 

GC 0053 

13 24 45.754 

-43 02 24.50 

19.43 

19.57±0.01 

26400 

2.22 

GC 0054 

13 24 46.435 

-43 04 11.60 

18.64 

18.84±0.01 

26400 

4.65 

GC 0058 

13 24 47.369 

-42 57 51.19 

19.15 

19.41±0.01 

24000 

2.69 

GC 0064 

13 24 50.072 

-43 07 36.23 

20.03 

20.11 ±0.02 

26400 

1.49 

GC 0065 

13 24 50.457 

-42 59 48.98 

19.68 

19.21±0.01 

26400 

2.17 


NOTE. — Summary of the new observations. Cluster identifications are listed in the first column, fol¬ 
lowed by the J2000 coordinates, apparent .R-band magnitudes used for target acquisitions, de-reddened 
apparent V- or r-band magnitudes (see §3.6), total integration times, and signal-to-noise ratios (S/N, 
see Sect. 2.2 for a definition). Table 1 is published in its entirety in the electronic edition of the Astro- 
physical Journal. A portion is shown here for guidance regarding its form and content. 


moil. Having performed these steps, the wavelength solutions 
were confirmed to be of high-quality by visually checking that 
they were smooth, as well as via the radial velocity errors in¬ 
ternal to the re-made slit geometry tables which showed val¬ 
ues of RVERR ~ 0.003 km s -1 . We note that these values 
are meant to confirm the accuracy of giwavecalibration and 
do not reflect our final, measured radial velocity uncertainties 
(see § 3). 

Using the final 11 sets of calibration data products, the 
recipe giscience provided the final, fully calibrated science 
frames from which individual ID spectra were extracted. Cus¬ 
tom Python scripts were used to clean the spectra of nu¬ 
merous residual cosmetic defects and to subtract the sky con¬ 
tribution from the spectra. To clean the spectra of cosmetic 
defects surviving the basic data reduction steps, the spectra 
were subjected to a median filtering algorithm and robust na- 
clipping. Each of the extracted spectra were visually inspected 
and the parameters of the median /no filters were adjusted to 
remove any significant detector cosmetics, while preserving 
the finer details of the spectra. Typically a median filter of 
gate size of 75 pixels followed by clipping points outside of 
4.5cr was sufficient to remove defects. 

2.3. Sky Subtraction 

To account for the sky contribution to each spectrum, we 
used the 13 GIRAFFE fibres dedicated to monitoring the 
sky contamination. These fibres facilitated uniform sampling 
across the field of view (see Figure 1). For each target, the sky 
contribution was taken to be the average of the three nearest 
sky fibres, inversely weighted by distance, thereby ensuring 
that only the sky nearest to each target was considered. The 
final sky spectra were determined individually for each of the 
117 targets and 11 OBs before being directly subtracted from 
each of the 1287 individual target spectra. Only then were the 
reduced, cleaned, and sky-subtracted spectra co-added to pro¬ 
duce the final data set. 

2.4. Data Quality Assessment 

The signal-to-noise ratios (S/N) for the final spectra were 
calculated considering the main spectral features used to esti¬ 
mate the line-of-sight velocity dispersions (see §3.1). Specifi¬ 
cally, these features are the Mg 6 and Fe 5270 Tick indices cen¬ 
tred at laboratory wavelengths of 5176.375 A and 5265.650 A, 
respectively (Burstein et al. 1984; Worthey 1994; Worthey & 


Ottaviani 1997). The S/N listed in Table 1 were calculated as. 


S/N 



( 1 ) 


where s,; and cr, are the mean and standard deviation of the 
flux over the continuum regions bracing the Mg/; and Fe 5270 
features as defined by 5142.625 A< A|)}™ corlt < 5161.375 A, 
5191.375 A < Aj5g 6|Cont < 5206.375 A, 5233.150 A < 
/ ^Fe 527 o,cont < 5248.150 A and 5285.650 A < A£ d 5270>cont < 
5318.150 A. Before calculating the S/N, each of the contin¬ 
uum definitions were shifted from the laboratory values to ac¬ 
count for known GC radial velocities, v r , or if unknown, they 
were shifted a posteriori according to our own v r measure¬ 
ments (see § 3.1). 


3. ANALYSIS 

3.1. Penalized Pixel Fitting 

Our line-of-sight velocity dispersion (EOS VD; a) measure¬ 
ments were carried out using the penalized pixel fitting (ppxf) 
code (Cappellari & Emsellem 2004). This code parametrically 
recovers the EOS VD of the stars composing a given cluster or 
galaxy spectrum by expanding the EOS VD profile as a Gauss- 
Hermite series. Using reasonable initial guesses for the ra¬ 
dial velocity (v r ) and cr, the best fitting u r>ppx f, cr ppx f, and 
Hermite moments / 13 , and /14 were recovered by fitting the 
cluster/galaxy spectrum to a library of template stars which 
had its spectral resolution adjusted to that of the FFAMES 
spectra. The fitting of optimal template spectra along with 
the kinematics serves to limit the impact of template mis¬ 
matches. An important feature of the ppxf routine is that dur¬ 
ing an iterative process, a penalty function derived from the 
integrated square deviation of the line profile from the best fit¬ 
ting Gaussian is used to minimize the variance of the fit. This 
feature allows the code to recover the higher order details in 
high S/N spectra, but biases the solution towards a Gaussian 
when S/N is low, as is the case for several objects in our sam¬ 
ple. For more details on the ppxf code, we refer to Cappellari 
& Emsellem (2004) 6 . 

Where possible, the input estimates for v r were quoted from 
Woodley et al. (2010a) or Woodley et al. (2007) which are 

6 ppxf and the corresponding documentation can be found at: http://www- 
astro.physics.ox.ac.uk/~mxc/idl/ 
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listed in Table 2. For GCs with unavailable v r , we used the 
IRAF 7 task rvcorrect to account for heliocentric velocity cor¬ 
rections, and fxcor to estimate v r . These estimates were used 
as our initial guesses for ppxf and do not need to be perfectly 
accurate since the cr ppx f measurements have no significant 
sensitivity to v r , as long as the initial guess is accurate to 
within a few tens of km s _1 (see also Taylor et al. 2010). The 
more refined iv, ppx f values were then adopted as our final esti¬ 
mates as listed in Table 2 and used for our S/N measurements. 

To account for any sensitivity that ppxf may have to the 
initial guess for a which may, for example, result in spurious 
solutions, we varied the input cr between 5 and 100 km s _1 in 
steps of 1 km s _1 for each spectrum. For each fit we used the 
entire 5152 < A/A < 5 365 range, including the Mg 6 and 
Fe features and all less prominent lines that are present. All 
GIRAFFE spectra have a velocity scale of 2.87 kms -1 , and 
we found an additive 4 th degree polynomial adequate for the 
purpose of estimating the continua. After varying our initial 
cr guesses, the adopted v r , ppx f and <r ppx f values correspond to 
the mean of the output sets of u r , ppx f and cr ppx f after being 
Kcr-clipped to remove outliers. The ppxf errors were taken as 
the mean of the output ppxf errors added in quadrature to one 
standard deviation of the Ktr-clipped results. While it is too 
cumbersome to present all of the ppxf fits in the present work, 
we direct the reader to the online-only appendix 8 to access 
the full suite of spectral fits for all targets in this sample and 
present a representative sample in Figure 2. 

3.2. Template Library 

The ppxf code measurements rely on a library of tem¬ 
plate stellar spectra capable of accurately replicating an inte¬ 
grated light spectrum when used in combination and Doppler 
shifted to account for velocity gradients within a stellar pop¬ 
ulation. We used a library of high-resolution synthetic spectra 
from the PHOENIX 9 collaboration (Husser et al. 2013). A 
synthetic spectral library was chosen over observed templates 
because current high resolution observed spectral libraries 
do not cover the wide range of stellar parameters that com¬ 
prise complex stellar populations. We therefore used a li¬ 
brary of 1 100 PHOENIX spectra covering the stellar param¬ 
eter ranges: 0.5 < log g < 4.0, 3 800 < T e g < 6 000 K, 
—4.0 < [Fe/H] < +1.0 dex, and —0.2 < [a/Fe] < +1.2 dex. 
For [Fe/H] = -4.0 and [Fe/H] =+1.0. The PHOENIX spec¬ 
tra are available only for [a/Fe] = 0.0 dex; however, we do 
not expect that this limitation affects our results significantly. 

3.3. Comparison with Previous Results 

The accuracy of ppxf when applied to the restricted ~ 200 A 
GIRAFFE wavelength coverage was tested by using the GC 
spectra of Taylor et al. (2010) as a comparison sample since 
they also derived LOSVD estimates in the same manner, but 
with much wider spectral coverage. To compare directly, the 
Taylor et al. spectra were constrained to the GIRAFFE wave¬ 
length range, and <7 ppx f estimates were obtained as described 
in the following. The input values for v r and a were fixed 
at those determined by Taylor et al., and the input a were 
varied around the known values by ± 10 kms -1 in steps of 1 

7 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron¬ 
omy, Inc., under cooperative agreement with the National Science Founda¬ 
tion. 

8 http://cdsarc.u-strasbg.fr/viz-bin/Cat/ 

9 http://phoenix.astro.physik.uni-goettingen.de/ 





Fig. 2.— Examples of the ppxf spectral fits (red lines) from a low-quality 
(top) to a high-quality spectrum (bottom). The reduced spectra (black lines) 
have been smoothed with a box size of 8 pixels, representative of the instru¬ 
mental FWHM of GIRAFFE. The inset blue plots show the residuals, beside 
which we list some output information for each target. From top to bot¬ 
tom this information is the GC identification, our adopted cr ppx f, iv.ppxf. 
reduced \ 2 , and the number of template spectra used in their calculation. 

km s -1 . The output values for a and v r are averaged and plot¬ 
ted as a function of input values in Figure 3 for comparison. 

The top panel of Figure 3 shows that the agreement for 
c ppx f is generally good, with most of the GCs clustering 
around the unity relation. The results of one GC (GC 0382) 
are not shown, as our new LOSVD value of 178.77+7.58 
kms -1 is unlikely to be reliable considering that it implies 
a dynamical mass of ~ 10 9 Mq within a half-light radius 
of ~ 2pc. Given that GC0382 has three independent v r 
measurements confirming it to be a member of NGC 5128 
(Woodley et al. 2010a, and the present work), and thus not 
a fore/background source, we adopt Taylor et al.’s cr ppx f = 
14.3 ± 3.2 km s -1 for the rest of the analysis. The few other 
outliers in Figure 3 correspond to the former study’s most un¬ 
certain GCs, so we prefer our er ppx f estimates since our tem¬ 
plate library has a significantly wider range of stellar param- 
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Fig. 3.— Comparison of the ppxf accuracy using the GC sample of Taylor 
et al. (2010) when the spectral range is restricted to that of the GIRAFFE 
spectra used in this work. The upper panel shows the a comparison and the 
lower panel shows the same for v r ■ In both panels, the Taylor et al. (2010) 
results are shown along the x-axes and the new results are shown along the 
y-axes. Unity relations are indicated as dashed lines. 


eters and much higher S/N ratio over the wavelength range 
used to estimate iv iPpx f and er ppx f. Meanwhile, the bottom 
panel shows that the agreement in rv !ppx f is excellent, with 
the scatter around the unity relation being consistent with 
the measurement uncertainties. The significant outlier corre¬ 
sponds to GC 0382, which we consider to be unreliable and 
defer to any previously derived v r estimates in the literature. 

Our new tv, PP xf estimates are listed in column four of Ta¬ 
ble 2, compared to literature values listed in column three. We 
note that GC0218, GC0219, GC0228 have v r measured for 
the first time (528 ± 1.90 kms -1 , 661 ± 2.70kms -1 and 
478 ± 19.40 km s -1 , respectively), all consistent with the 
541 kms -1 systemic velocity of NGC5128 (Woodley et al. 
2010a). There were 15 GCs for which ppxf was unable to pro¬ 
vide v r estimates, including GC0261 and GC0315, leaving 
them still as v r unconfirmed members of NGC5128. Thus, 
Table 2 lists new accurate v r estimates for 125 NGC5128 
GCs, including three first-time measurements. 

We compare the new tv iPpx f estimates to literature radial 
velocities (rv,iit) in Figure 4 where v r ^t and w r . ppx f are shown 
along the x- and y-axes, respectively. The Woodley et al. 
(2007) and Woodley et al. (2010a) catalogues provide the 
most comprehensive collections of NGC5128 GC radial ve¬ 
locities to date. For the comparison we adopt the weighted- 
average values listed in Woodley et al. (2010a) where pos¬ 
sible. If there exists only a single Woodley et al. (2010a) 
value measured from an individual spectmm, then we adopt 
the Woodley et al. (2007) estimates, unless the former agrees 
significantly better with our new measurements. Figure 4 
shows generally good agreement within the literature uncer¬ 
tainties, with a single notable exception being GC0095. For 
this GC, we prefer the literature measurement of u r j it = 
374 ± 34kms -1 over our cr ppx f = 826 ± 7.8 kms -1 because 
visual inspection of the corresponding Doppler shifts shows 


750- 


650- 


03 5 5 0- 

a 

450- 


g;350- 


250- 


150- 



50- 


50 150 250 350 450 550 650 750 850 

- 1 - 


<V,iit t km s ] 


Fig. 4.— Comparison between the radial velocities newly estimated by 
ppxf, ^r,ppxf t0 those previously measured from the literature, The 

comparison is generally good, as indicated by the dashed line showing the 
unity relation. However, our measurements are in most cases of significantly 
higher quality, as illustrated by the u rjPpx f error bars, which for our dataset 
are often on the order of or smaller than the symbol size. 


better agreement with the laboratory wavelengths of the spec¬ 
tral absorption features when using the former value. Given 
this discrepancy, and the resulting uncertainty of the derived 
(jppxf, we drop this object from the analysis. In any other case 
where v r are discrepant we prefer our ry iPpx f due to smaller 
uncertainties. 

Despite obtaining reliable v r estimates for almost all of our 
targets, there were 25 GCs for which we consider cr ppx f to 
be unreliable either due to uncomfortably large error bars, 
or simply a failure to derive an estimate at all. We therefore 
drop these targets from the subsequent analysis and carry on 
with the remaining 115 new er ppx f estimates, including the re¬ 
analyzed Taylor et al. GCs. 

3.4. Structural Parameters 

We took 2D projected half-light radii, r^, and concentration 
parameters, c, from the three sources as listed in Table 3. The 
majority of the values listed for rh and c are, where available, 
from Iordan et al. (2015, in prep.), based on HST data, and 
are otherwise taken from Gomez et al. (2015, in prep.) based 
on IMACS data taken under exceptional seeing (< 0.5"). De¬ 
spite the sub-arcsecond seeing conditions, the marginally re¬ 
solved nature of most GCs did not allow for accurate c esti¬ 
mates by Gomez et al., so for many, a typical value of 1.48 
was assigned. Moreover, three extremely small IMA C.S'-based 
j'i, estimates exist where there is no 11ST imaging available; 
however, these GCs (GC0085, 0333, 0429) are not resolved 
in the images and thus we drop them from the analysis and 
continue with the remaining 112 CSSs. The re-analyzed clus¬ 
ters of Taylor et al. (2010) use the same parameters as in that 
paper, namely those derived in Harris et al. (2002a). For the 
latter GCs, the 77 , errors listed are adopted from the same pa¬ 
per. At the time of writing, no errors were available for the 
Gomez et al. and Iordan et al. sizes, so for these clusters we 
adopt values of 0.43 pc, which is the average size error of the 
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Taylor et al. (2010) sample, and while representative for the 
new HST measurements, may underestimate the IMACS er¬ 
rors. Following Harris et al. (2002a), we assume 0.15 for the 
errors on all c, noting that these probably underestimate the 
true IMA C.S'-hased measurement errors as well. In any case, 
although assigning error estimates in this manner is not opti¬ 
mal, we do not expect it to affect our main results significantly 
given the dominance of u on the uncertainties of our dynami¬ 
cal mass estimates (see § 3.6). 

3.5. Aperture Corrections 

The < 7 pp X f values listed in Table 2, while generally accu¬ 
rate, are not appropriate to use when estimating dynamical 
masses. Several effects, both observational (seeing, target dis¬ 
tance, etc.) and instrumental (spectral/spatial resolution, sam¬ 
pling, etc.), may conspire to affect how representative the 
light entering a given fibre aperture may be of objects simi¬ 
lar to massive GCs and UCDs (Mieske et al. 2008a). In our 
case the 1.2" diameters of the FLAMES fibres correspond to 
~ 22 pc at the distance of NGC 5128, so contributions from 
stars outside of the core region may skew the <T ppx f measure¬ 
ments to lower values compared to a estimates corresponding 
to smaller radii. Here we describe our approach to correct our 
measured a ppx f estimates to values representing both the GC 
core regions (oq) and a values within the GC half-light radius 
(f J l/ 2 )- 

We used the cluster modeling code of Hilker et al. (2007), 
described in detail in Mieske et al. (2008a), to correct for 
any aperture effects and determine estimates of <ro and <Ji /2 
based on our measured rr ppx f. This code uses the basic struc¬ 
tural data (in our case />, and c) that defines a cluster’s light- 
profile to generate a 3D King (1966) stellar density profile 
from which an N-body representation of the cluster is created 
in 6 D (position, velocity) space. Each simulated particle is 
convolved with a Gaussian corresponding to the true seeing 
FWHM (see § 2) and a light profile is generated from which 
the velocity dispersion profile can be obtained. 

Using this code, we modeled each of our clusters with 10 5 
particles and binned them radially in groups of 10 3 . The 3D 
velocity information of each subgroup was used to derive a 
profiles according to (^) 3D = 3cr" where ( v )!d is the square 
of the mean of the 3D velocities. To account for the inherent 
stochasticity of the modeling, the median of the inner-most 
five subgroups, or 5% of the modeled stellar population, was 
adopted as cto, while all particles inwards of r/, were used to 
calculate <ri / 2 for each GC. This process was repeated three 
times per GC. The first set of models used the measured <7 ppx f, 
Th, and c as inputs to provide our adopted <r 0 and cr i/ 2 > while 
for the other two iterations we added or subtracted the errors 
for the three quantities in order to maximize or minimize the 
modeled cr estimates, respectively. We then adopted the dif¬ 
ferences between the upper/lower bounds and the output cr 
values as the corresponding errors. Table 2 lists the resulting 
cr 0 and estimates including our uncertainties (columns 8 
and 9), alongside the model velocity dispersions correspond¬ 
ing to the FLAMES apertures (cr mj0p ) and the global values 
The accuracy of the code is verified by the very good 
agreement between the predicted cr ppx f and measured cr m op 
at the fibre aperture size. 

3.6. Star Cluster Masses and Mass-to-Light Ratios 

One of the most direct methods to estimate the dynami¬ 
cal mass (Addyn) of a single-component compact stellar sys¬ 
tem is by the use of the scalar virial theorem (e.g. Binney 


& Tremaine 1987) of the form originally derived by Spitzer 
(1969), 

<2) 

if one assumes a dynamically relaxed cluster, sphericity, and 
isotropic stellar orbits. While this is among the most com¬ 
monly used dynamical mass estimators, it has been shown that 
the “half-mass” (Ad 1 / 2 ) or in other words the dynamical mass 
corresponding to that contained within the 2D projected half- 
light radius is more robust against stellar velocity dispersion 
anisotropy. This feature makes Adi / 2 an overall more robust 
mass estimator for dispersion supported systems. We estimate 
Ai\/ 2 via the form derived by Wolf et al. (2010), 

M I/ , = 4<%tt=,93o(A^L)(^)M 9 (3) 

where <ji os is the luminosity-weighted LOSVD, in our case 
aperture corrected to g \/ 2 . 

Applying Equation 3 to all the GCs with available 7\ / 2 , 
r> ( , and c provides Ai\/ 2 estimates for a total of 112 of 
NGC 5128 star clusters, 89 of which are first-time measure¬ 
ments, in particular at faint absolute luminosities (see Sec¬ 
tion 4.1 and Figure 7). We find in our star cluster sample 
Adi / 2 estimates ranging from the low-mass end. Ad 1 / 2 ,min = 
3. 7^37 • 1O 4 M 0 (GC0031), to the highest-mass object, 
GC 0365, with Ad 1/2 , max = 7.41^;^ • 10 6 M 0 , with a sam¬ 
ple median Adi / 2 of 3.47-lO 5 M 0 . By assuming that mass 
follows light, these masses translate into total mass estimates 
of Adtot.min = 7 .iAlAf 1 • 10 4 M 0 representative of the lower 
range of GC masses. Ad to t,max = l-48to'gg • 10' M 0 consis¬ 
tent with UCD masses, and a median Ad tot = 6.94- 10 5 Mq. 

Alternatively, we use uo to estimate Add yn via Equation 2. 
By doing this, we lose the benefit of isotropy independence, 
but gain freedom from the underlying assumption that mass 
follows light. Encouragingly, we find very similar masses us¬ 
ing either estimator at the low-mass tail of the GC mass distri¬ 
bution, with Addyn,min = 7.6ly g 7 • 10 4 Mq, and at the high¬ 
est mass with Add y n,max = 1 - 53+0 15 ‘ 1 O 7 M 0 and a me¬ 
dian Add y n of 7.04 • 10 5 Mq. Interestingly, for total masses 
higher than ~ 1.0 • 10 6 M 0 we find that the relative difference 
between the estimators becomes more significant for certain 
clusters. 

Figure 5 shows the relative difference between the total dy¬ 
namical masses derived using 2 • Adi / 2 and Ad d yn as a func¬ 
tion of Addyn for each GC in our sample. Add yn essentially 
predicts similar (within ~ 1 — 5%) masses for “normal” GCs, 
(< 10 6 M 0 ; typical of Local Group GCs). Above this thresh¬ 
old the discrepancy between the two mass estimates becomes 
more pronounced for certain clusters, reaching up to ~ 18% 
higher Add yn - Altogether, this comparison suggests that for 
typical GC masses, i.e. in the range ~ 10 5 —10 6 M 0 , where 
cluster masses are completely dominated by baryonic mate¬ 
rial, Adi /2 is the more robust measure of the total GC mass. 
Above this mass range, and outside of r> ( any kinematical trac¬ 
ers arising from, e.g. non-equilibrium configurations or/and 
dark gravitating mass components, may introduce biases in 
the mass estimates. 

We calculate the dynamical mass-to-light ratios evaluated 

1 /2 

within the half-light radius (T v - ) by dividing 2-Adi / 2 by the 
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TABLE 2 

Star Cluster Kinematics 


ID 

R g c 

Vr, lit 

4V,ppxf 

Oppxf 

Om,ap 

&m,gl 

&0 

c 1/2 

Ref. 


[arcmin] 

[km s -1 ] 

[km s -1 ] 

[km s -1 ] 

[km s -1 ] 

[km s _1 ] 

[km s -1 ] 

[km s _1 ] 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

GC 0028 

11.30 

558±97 

512±1.80 

8.95T2.10 

8.97±|‘?f 

8.60ii;®° 2 

Q 07+2.26 

9-o7_2 25 

9 32"*" 2 ’ 27 

^• JZ -2.25 

1 

GC 0031 

10.63 

595T202 

552±2.80 

3.69T4.35 

3.«tS;i5 

3 34+4.05 
-3.53 

4-12ll.II 

3 94+ 4 68 
-3.93 

1 

GC 0048 

11.36 

509±15 

481±1.10 

7.61±1.40 

7-62li; 3 3 8 8 

7-08±Ll! 

8-25±l;®| 

7-9311S 

1 

GC 0050 

B.03 

718±16 

732±0.70 

10.09±0.90 

10 .01±g;®| 

9.oi±g;5i 

10 65" 1 " 0 ' 82 
iu.oj _ 0 96 

10 .081“;“ 

1 

GC 0052 

7.89 

276±59 

199± 15.25 

1.09±51.85 





2 

GC 0053 

7.75 

503±17 

505T1.00 

10.50±1.20 

10.55±i;|g 

9.951?;® 3 

nmtlH 

11 31+ 1 - 29 

1 

GC 0054 

8.12 

736±58 

734±1.50 

9.31T1.70 

9 30 +1 ' 72 
1.68 

8.16±LI! 

9-60 till 

9 23“*" 1,63 

3 

GC 0058 

8.06 

685±43 

684±0.80 

7.60=bl.l0 

7 61" 1 " 1 ' 06 

7 47+0-92 
//f/ -1.03 

8.50±LJ| 

8-24ll; 2 3 ? 

2 

GC 0064 

9.42 

594±36 

574±2.90 

16.29=b2.90 

16 27 -1 " 2,89 
io.z / _ 2 89 

15 29 +2,4 ° 
lJ -* y - 2.50 

17 79+ 3 35 

1 / - ,y -2.95 

17.0811:5* 

1 

GC 0065 

6.92 

331±47 

278±1.30 

9.78=bl.62 

9 78“*" 1,58 
y ‘ /0 —1.64 

o | o + 1.26 

yLD -1.37 

11 25" 1 " 1 ’ 76 

10.55lL.?| 

2 


NOTE. — Kinematical data for the NGC 5128 star clusters. Cols. 1 and 2 list the cluster IDs and projected galacto-centric radii respectively, cols. 3 
and 4 list radial velocities, and cols. 5-9 list <r measured with ppxf and values which have been aperture-corrected to various cluster radii (see § 3.5 
for details). Where available, all v r ,ut values are taken from Woodley et al. (2010a) corresponding to their “mean” values otherwise we adopt the best 
matches from either Woodley et al. (2007) or the estimates from Woodley et al. (2010a) that are based on individual spectra. Table 2 is published in its 
entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. 

References for i> r ,iit- (1) Woodley et al. (2010a); (2) Woodley et al. (2007), their mean; (3) Woodley et al. (2007), LDSS2; (4) Woodley et al. (2007), 
VIMOS; (5) Woodley et al. (2007), Hydra. 


TABLE 3 

Star Cluster Structural Parameters 


ID 

•Mi /2 

A4dyn 

Tl /2 

T dyn 

rh 

c 

Ref. 


[M 0 ] 

[M e ] 

[M©/L q ] 

[Mq/Lq] 

[pc] 



( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

GC0028 

0 25 - *" 0,12 
u - zd _ 0.12 

°- 5 °la24 

3 42+ 1,71 

j.^ z _l 70 

3 45+1-66 

2.96=b0.43 

1.48 

1 

GC0031 

0 Q4+° 09 
u - u ^_0.04 

O-OSlo'.oi 

0 49 + 1 16 
u/fy -0.48 

0 50 +1 16 

0.30-0.49 

2.58=b0.43 

1.48 

1 

GC0048 

0 24 - *" 0 09 
U - Z ^-0.09 

0-48±g;ig 

2 3Q+0-98 

zz,u -0.96 

2 54+0-99 
z - J ^-0.98 

4.08=b0.43 

1.65 

2 

GC0050 

0 66”^°' 13 
u - oo - 0.12 

l-39lo'.26 

3 43 +O .68 

•3-^-0.63 

3 5Q+0-62 
3.3V-0.70 

7.03T0.43 

1.56 

2 

GC0052 





1.02=b0.43 

1.48 

1 

GC0053 

a 0-3+0.09 

u - DD -o.io 

0 . 68 l°;i* 

3 ft] +102 

J-Oi_1.08 

3 76+ 1 ' 00 
3./°_i 09 

2.74±0.43 

1.83 

2 

GC0054 

0 71+9-25 

u - /a -0.27 

1 43+0-56 

i-^J-0.56 

4 04+ 1 -45 
^• u ^-1.56 

4 iq+ 15 9 

9.01T0.43 

1.65 

2 

GC0058 

0 12 - *" 0,95 
u - iZ -0.05 

0.241°;°® 

1 14+0-43 
1 - i ^-0.45 

ii4l8;lt 

1.89=b0.43 

1.43 

1 

GC0064 

l- 02 t°jl 

2.08l°;*i 

18.651^1 

18.97 t 7 6 f 2 

3.77=b0.43 

1.57 

1 

GC0065 

0-28lo.il 

°- 6 °la23 

7 7 ^ +0.81 
z - z:> -0.84 

7 to+0.84 
z -3O_ 0 .92 

2.71 ±0.43 

1.99 

1 


NOTE. — Structural parameters and dynamical masses of our GC sample. GC0082, GC0107, GC0214, 
GC 0219, GC 0236, GC0262, GC0274, GC0417, GC0418, GC0420, GC0439, GC0435, and GC0437 are 
based on de-reddened r' band data. Table 3 is published in its entirety in the electronic edition of the Astro- 
physical Journal. A portion is shown here for guidance regarding its form and content. 

References for rn and c: (1) Gomez et al. (2015, in prep.); (2) Jordan et al. (2015, in prep.); (3) Harris et al. 
(2002a). 


total, de-reddened V-band luminosity, calculated as, 

Ijy = iQ-0.4(V o -(m-M)°-Mv,a) ( 4 ) 

where My ;0 =4.83 mag. We also calculate dynamical mass- 
to-light ratios based on <r 0 (T^J n ) by dividing A4d yn by 
Ly■ Most of our sample have apparent L-band magnitudes 
provided in the Woodley et al. (2007) catalogue (see also ref¬ 
erences therein), for which we list the de-reddened values in 
Table 1. We account for foreground reddening on an individ¬ 
ual basis, with no attempt to correct for extinction internal to 
NGC 5128, by using the Galactic Extinction and Reddening 
Calculator 10 with the galactic reddening maps of Schlafly & 
Finkbeiner (2011). Where i '-band magnitudes are not avail¬ 
able, we list de-reddened r'-band magnitudes from Sinnott 
et al. (2010). If no photometry in the previously mentioned 

10 http://ned.ipac.caltech.edu/forms/calculator.html 


filters is available, we base our T y estimates on the /?-band 
magnitudes from the acquisition images, assuming in all cases 
a conservative photometric error of 0.1 mag. 

1 /2 

We find a large spread in the corresponding T y ranging 
between 0.49 M 0 Ly ^ for GC 0031 up to 64.47 M 0 i+^ for 

GC0225, with a sample median T(/ 2 of 3.33 MqL^q. Sim¬ 
ilarly, we find 0.50 < Yy yn < 66.61 M 0 I+q with a slightly 
higher median value of 3.44 M 0 I/y 0 . While these results 
are smaller by 0.47 and 0.36 M 0 Ly 0 than the median of 
3.8 M Q Ly) Q found by Taylor et al. (2010) for NGC 5128 
GCs, they are notably higher than the median Yydyn of 
2.2 ± 0.3 for Milky Way (MW) GCs (McLaughlin 2000; 
McLaughlin & Fall 2008) and Tu.dyn = 1-37 ± 0.28 for 
M31 GCs (Strader et al. 2011). This result should perhaps 
not be too surprising as our sample is biased toward GCs at 
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the bright end of the GC luminosity function (GCLF). Thus, 
we are most likely not including many GCs with typical Lo¬ 
cal Group GC masses so that our sample is biased to GCs 
above the ~ 2 • 10 6 Mq threshold where Yy/ n begins to 
rise dramatically (e.g. Hagegan et al. 2005; Kissler-Patig et 
al. 2006; Mieske et al. 2006, 2008a; Taylor et al. 2010). Con¬ 
versely, our sample includes not only those of Taylor et al., 
but many fainter GCs, thus reaching well below the afore¬ 
mentioned threshold and biasing our medians toward slightly 
lower values compared to previous studies. 

4. DISCUSSION 

Correlations between Yp vn , Add yn , and absolute magni¬ 
tude (My) can provide important information on the dynamic 
configuration and baryonic makeup of star clusters. To inves¬ 
tigate these relations, Ty yn is shown as a function of Addyn 
and My for each of our sample GCs in Figures 6 and 7 re¬ 
spectively, with the color shading parametrizing GC i'h- 

Given that the following discussion hinges strongly on the 
features seen in Figures 6 and 7 having astrophysical explana¬ 
tions, we first considered several systematic effects and per¬ 
formed corresponding tests to check whether these could bias 
our measurements and artificially generate the observed re¬ 
sults. The description of these tests, including detailed checks 
for data analysis biases, target confusion, correlations with 
galactocentric radius (R gc ), and/or insufficient background 
light subtraction is provided in the Appendix. In summary, 
none of the tested effects are likely to explain the observed 
features in Figure 6 and 7, and thus we consider astrophysical 
explanations in what follows. 

4.1. Addyn vi. T^ vn Relations 
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Fig. 5. — Relative difference between the traditionally used dynamical 
mass Mdyn and twice the half-mass, i.e. -Vl,i yll — 2 -.M qq, as a function of 
dynamical mass, .Mdyn- If mass follows light, then the total GC mass should 
be approximated by twice A4 i/2 > hence, the difference being close to zero 
(dashed line). Here, Mdyn predicts somewhat higher masses than 2- M 1/2 
for GCs with A4dyn A10 6 M 0 , indicating that for these clusters, mass may 
not strictly follow light and either dark gravitating mass components or non- 
equilibrium configurations may become necessary to explain the discrepancy 
in dynamical mass estimates of these clusters. 


A number of interesting features shown by Figure 6 regard¬ 
ing the mass, size and mass-to-light ratios of our sample GCs 
are of note. We see a clear bifurcation in the Ty Vn -Add yn 
relations at Add yn ^ 10 6 Mq, with two well debited se¬ 
quences of GCs showing distinct positive slopes. GCs below 
Addyn ~ 10 6 Mq do not seem to follow either of these two 
relations, appearing to have mass-to-light ratios in the range 
of ~ 1—5 with no particularly well debited correlation. In gen¬ 
eral, the circular points in Figures 6 and 7 indicate a smooth 
transition from GCs with masses Add yn ~ 10 5 — 10 6 M 0 to 
those with > 10 6 Mq that follow the two sequences. 

To compare our measurements with similar data of Local 
Group GCs we overplot measurements from McLaughlin & 
van der Marel (2005) for Milky Way (MW), Large Magel¬ 
lanic Cloud (LMC), and Small Magellanic Cloud (SMC) GCs, 
along with data taken from Strader et al. (2011) for M31 
GCs. It is important to point out here that while the Strader 
et al. measurements are based on direct kinematical measure¬ 
ments, the same cannot be said for the McLaughlin & van der 
Marel (2005) data since they are based op estimates extrap¬ 
olated from modelled light prohles. Thus, while we use the 
McLaughlin & van der Marel data as a large, homogeneous 
comparison dataset, a true comparison cannot be made un¬ 
til Addyn measurements can be made for a large sample of 
MW GCs based directly on stellar kinematics. With that said, 
both samples align well with the bulk of our NGC 5128 GC 
sample but extend to signibcantly lower masses and fainter 
luminosities. From the comparison with the majority of the 
Local Group GC sample, we conclude that the NGC 5128 GC 
sub-sample with Add yn 10 6 M 0 and Tj y yn < 5 MqLq 1 can 
be regarded as “classical” GCs similar to those found in the 
Local Group that follow the well-known GC “fundamental 
plane” relations (Djorgovski 1995; McLaughlin 2000) where 
non-core-collapse GCs in the MW show almost constant core 
mass-to-light ratios of 0 = 1.45 MqLq 1 . 

It is difficult to establish whether the fundamental plane re¬ 
lations strictly hold for the “classical” NGC 5128 sub-sample 
as the core surface brightness values of these clusters are 
not accessible, even with HST imaging. Even so, it can be 
seen in the middle panel of Figure 8, which compares the 
mass/light properties of our NGC 5128 sample with the Lo¬ 
cal Group GCs that the three distributions show a very simi¬ 
lar rise up to Tj;’ rl ss 2 MqLq 1 ; above which the NGC 5128 
sample dominates. With a MW/SMC/LMC sample that under¬ 
represents the peak of the GCLF at Ly,© ~ 2 • 10 5 L 0 (see 
bottom panel of Figure 8), it is beyond the scope of this work 
to verify whether a larger sample of MW/SMC/LMC GCs 
would “hll in” the distribution shown by our NGC 5128 sam¬ 
ple. The M31 sample, on the other hand, tends toward lower 
Addyn, despite sampling a similar luminosity range as our 
NGC 5128 data. We thus, for now, consider the NGC 5128 
GC sub-sample shown by the circular points in Figures 6 and 
7 to simply represent the GCs that populate the Local Group. 

More generally, we note that most NGC 5128 GCs have 
Ty™ greater than the median value of MW/SMC/LMC GCs 
(T\f\y ~ 2.2 MqLq 1 ) and signibcantly higher values than 

for M31 GCs (Y^^ «1.37 MqLq 1 ), although the latter may 
be due to insufficient aperture corrections for extended GCs 
(see discussion in Strader et al. 2011). With a median Tj, " of 
3.44 MqLq 1 , it is tempting to suggest a fundamental differ¬ 
ence between Local Group GCs and those in NGC 5128. We 
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Fig. 6.— T'j 1 '" 1 as a function of .Vl,i yll . with the color shading indicating the half-light radius of the clusters. Three groups of GCs corresponding to “classical” 
GCs, and to the “low-” and “DSC” sequence GCs are indicated by circular, pentagonal, and square symbols, respectively (see legend). The empirical power- 
law relation of Taylor et al. (2010) is shown by the thin dashed black line, while steeper power-law fits based on the two Y?/ n branches are shown by the 
colored thick dashed lines. The steepest relation (red) is based only on the square points which have been isolated based on their faint absolute magnitudes 
(My > —8.5 mag) and high mass-to-light ratios (Tj, yn > 6.0 Mq Ly y ). Meanwhile, the blue dashed line shows the results of fitting only the pentagonal points 
with Ty° < 10.0 A7q L y 1 , and ,Vl,i yll > 2 ■ 10 6 Mg. The solid black line and dark gray shading illustrate the observational limit to our data, while the light 
gray shading indicates an approximate lower mass limit for selected GCs in our spectroscopic sample. Several interesting objects that are discussed in the text are 
labeled with their catalog numbers. We plot corresponding measurements for Local Group GCs as grey stars (MW/SMC/LMC; taken from McLaughlin & van 
der Marel 2005) and grey triangles (M31; taken from Strader et al. 2011). 

do not strictly support this notion, as a more straight-forward 
interpretation is that we are simply sampling part of the 
GCMF that is inaccessible in the Local Group due to a dearth 
of known GCs above 10 6 Mq. Indeed, the top panel of Fig¬ 
ure 8 demonstrates that the A4 dyn distributions of the Local 
Group and NGC5128 GC samples are very dissimilar; the 
highest values are reached by NGC5128 GCs, followed by 
M31 and the MW/SMC/LMC. Above M dyn « 5 • 10 5 M e , 
the NGC5128 sample is well represented up to ~ 10 7 Mq, 
while Local Group GCs are more populous toward the lower 
tail of the GCMF. Altogether, Figure 8 suggests that the 
NGC 5128/Local Group GCSs may have fundamentally dif¬ 


ferent GCMFs, given the similarly sampled GCLF (bottom 
panel), something that can be tested when similar NGC 5128 
data probing fainter magnitudes becomes available. 

Having addressed the main, “classical” body of GCs in Fig¬ 
ure 6, we now turn to the two distinct high-.Ad,i yri sequences. 
In the following we refer to GCs with T^ yn > 6 and luminosi¬ 
ties fainter than My = —8.5 mag (see Figure 7) as members 
of the “dark star cluster” (DSC) sequence (red dashed line in 
Figure 6) due to their potential connection to DSCs predicted 
by theory (see §4.3.2; Banerjee & Kroupa 2011). Those with 
A4 dyn > 2 • 10 6 Mq and Yy yn < 10 that follow the shal¬ 
lower T', l /"-A / (,iy I , relation (blue dashed line in Figure 6), we 
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Fig. 7.— Trf" as a function of My , with cluster half-light radius indicated by the color bar. Symbols are the same as in Fig. 6, with the same objects indicated. 
The vertical light-grey line and the shaded zone (My > —7.4 mag) indicate the photometric limit of our spectroscopic sample selection. 


refer to as members of the “dwarf-globular transition object” 
(DGTO) branch, since these objects encroach upon the struc¬ 
tural parameter space of the DGTOs reported by Ha§egan et 
al. (2005). 

Interestingly, the two objects omitted by these criteria, 
GC0217 and GC0239, lie intermediate between these two 
sequences and so we refer to them as “intermediate-Y 1 ] 1 / 11 ” 
star clusters (see Figure 6). Their properties may indicate ei¬ 
ther an evolutionary connection to one of the sequences, or 
perhaps represent a separate population that is simply not 
well sampled by our data. In terms of structural parameters, 
the mean of the DSC sequence ({ryjhigh = 5.10pc) is 
marginally smaller than that of the DGTO branch ((r/j)i ow = 
6.67pc), while also smaller in the median (4.01 and 5.59pc, 
respectively). Meanwhile, the mean galactocentric radii of 
the two populations are (i? gc ) high = 6.97' and (i? gc )iow = 
9.11'. Welch 2-sample tests yield that the mean differences 
in rh and R gc are not statistically significant, with p-values of 
0.27 and 0.15, respectively. 


4.2. Properties of the Y^ l : vn -Aldyn Sequences 

To probe the properties of the two Tj;’ n -,M,i y i l sequences, 
we fit empirical power-law relations of the form, 

Ty n a M2 yn (5) 

to approximate the data. In a similar analysis of dispersion 
supported CSSs, including a subsample of the GCs consid¬ 
ered here, Taylor et al. (2010) found a value of a = 0.24±0.02 
to fit their data, connecting “classical” GCs to more massive 
systems like UCDs and dwarf elliptical galaxies. This rela¬ 
tion is shown in Figure 6 by the thin dashed black line and is 
too shallow to fit either of the Ty yn -yVldyn sequences of the 
present GC sample. To better represent the data, we instead 
make an effort to fit power-laws to the sequences individu¬ 
ally and find that each is well approximated by distinct, tight 
relations as described in the following. 

The dashed blue line in Figure 6 shows an approximation to 
the DGTO sequence with a power-law slope of oi = 0.33±0.04 
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Fig. 8 .— Comparison between the Aidyn ( to P panel), Tj iy " (middle 
panel), and Ly (bottom panel) distributions for the NGC 5128 objects of 
the present work (red histograms), and Local Group GCs (MW/SMC/LMC; 
(blue histograms), and M31; (green histograms)). Opacity has been lowered 
to show where the distributions overlap. While the MW/SMC/LMC sample 
is too small to make strong comparisons to the other two distributions, it is 
clear that the NGC 5128 sample shows higher median Aidyn and Tr7 n than 
the M31 sample, despite probing a similar range in Ly. The total sample 
sizes are indicated in parentheses in the legends. 


(pentagons in all relevant Figures). This relation fits the data 
quite well from the high mass GCs down to ~ 10 5 Mq, which 
represents a value more typical of Local Group GCs. Mean¬ 
while, the DSC sequence (square points in all relevant Fig¬ 
ures), shown by the red dashed line in Figure 6 with a steeper 
slope (a = 0.79 ± 0.04), seems to be created by a funda¬ 
mentally different collection of objects. Interestingly, we find 
two LMC GCs (NGC 2257 and NGC 1754, see Figure 6) that 
appear to align well with the DSC sequence. While no strong 
statements can be made about only two objects, their exclu¬ 
sive presence around a currently interacting satellite of the 
MW marks an interesting starting point to investigate any con¬ 
nection to the DSC sequence. 

Ty yn is plotted as a function of My in Figure 7, which 
shows that the DGTO sequence is composed exclusively of 
the brightest GCs of the sample. Thus, the DGTO sequence 
may simply be explained by these GCs representing the high- 
luminosity tail of the GCLF. On the other hand, GCs on the 
DSC sequence are fainter than DGTO GCs by A My «1 mag, 
making them similar in luminosity to the average GCLF turn¬ 
over magnitude found in many GC systems. Furthermore, it 
can be seen that the range of /•/, for the objects on the DSC 
sequence is not dramatically different from that of the “clas¬ 
sical” or DGTO GCs. Collectively, the similarity shown by 
these DSC objects in luminosity and size to other GCs in 
many GCSs likely explains why they have not been identified 
before in other galaxies, as they are only remarkable in their 
stellar dynamics properties. Regardless, as this is the first time 
a clear distinction between two such groups of CSSs has been 
made, this naturally leads to the question of whether these ob¬ 
jects should be called GCs at all. 


4.3. Possible Origins of the Sequences 

Having shown artificial biases to be unlikely drivers of our 
results (see Appendix), the following discusses several astro- 
physical mechanisms that may generate our observations. 

If the DGTO sequence is made up of the brightest “clas¬ 
sical” GCs, the Tj 1 -'’ 11 ~ 5 — 9 values shown by some are 

still perplexing in that GCs with T^ n > 5 require addi¬ 
tional explanations beyond being the extension of the “clas¬ 
sical” GCLF. Effects that can mimic higher than usual 
include non-equilibrium dynamical processes (e.g. rotation, 
pre-relaxation, young stellar populations, tidal disruption) or 
an exotic IMF. We find several DGTO sequence GCs with el- 
lipticities, e > 0.25, indicative of a non-equilibrium dynami¬ 
cal state, such as rotation. Other GCs show young (< 8 Gyr) 
ages (Woodley et al. 2010b), and thus may not be fully re¬ 
laxed. Additionally, for four objects on this DGTO sequence 
(GC 0041, GC 0330, GC0365, and GC0378) Harris et al. 
(2002a) found evidence for extra-tidal light contributing to 
their surface-brightness profiles in excess of their King model 
fits. For the remainder, the possibility of a particularly bottom- 
heavy IMF (e.g. Dabringhausen et al. 2008; Mieske et al. 
2008b) could explain their elevated Ty yn estimates. 

4.3.1. Globular Cluster Rotation 


To probe how non-equilibrium states could explain the el¬ 
evated Ty yn , we investigate the possible impact of rotation 

on our mass estimates. Treating the observed Yj lyn as be¬ 
ing strictly the result of the GCs exhibiting a mass excess at 
a given luminosity with respect to the median T^w = 2.2, 
we calculate for each of our sample GCs the amount of “ex¬ 
tra” mass within r h , or AM 1/2 = ^Xy yn -T^) • L f- = 

^T^ iyn —2.2^j • This mass component then needs to be 
accounted for by the effects discussed above to explain the 
elevated Yj*. yn values. 

Making the naive assumption that AM i / 2 is entirely due to 
rotation, then in the general situation where the rotation axis is 
aligned at any angle i with respect to the observational plane, 
stars at ry, would require circular velocities of at least. 


v sm i ~ 


^ AM 1 / 2 G y /2 


(6) 


To investigate how rotation could explain the Yj lyl esti¬ 
mates, the left panel of Figure 9 shows a comparison between 
v sin i and the average stellar velocity in GCs computed from 
the random stellar motions via ct\ /■> V-i. The symbols are the 
same as in Figures 6 and 7, but with “classical” GCs sim¬ 
ply shown as dots. The dashed line in Figure 9 indicates the 
boundary at which rotation is equal to the random stellar mo¬ 
tion component required to be consistent with AM j /2 . It must 
be acknowledged that the errors bars shown on Figure 9 make 
it impossible to definitively discuss the dynamical configu¬ 
rations of our sample. On the other hand, the lack of co¬ 
mingling between the different groups suggests that the fea¬ 
tures seen are probably not solely due to systematic errors, 
thus they can still be used to make general statements about 
the populations as wholes, and we proceed with that in mind. 

A GC that has a rotation to random motion component ra¬ 
tio, v sini/cr lj / 2 y/3 > 1, requires circular velocity speeds that 
would destabilize the system if only rotation is to explain its 
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Fig. 9.— (Left panel): Average measured orbital velocities of cluster stars shown as a function of the minimum circular velocities necessary for rotation alone 
to account for the observed mass-to-light ratios. The dashed line represents the critical velocity (assuming sin i = 1. or the maximum observable velocities) 
at which a) rotation alone is able to explain the elevated T^ yn of an object, or/and b) an object of a given mass becomes unstable against rotation, consistent 
with a dynamical non-equilibrium state. Symbol shapes are as in Figures 6 and 7, but with the “classical” GCs shown as points and the two intermediate- 
n GCs (GC0217 and GC 0239, see Section 4.1) shown as gray symbols. (Right panel): v sin y 2 as a function of the dynamical half-mass, .VI 

for the same objects. There are correlations for the objects on the DSC and DGTO sequences, with a weaker correlation shown by the “classical” GCs with 
v sin %I\/ 3(T | /2 > 0.2. indicating that GCs become more rotationally supported at higher masses. 


AM |/ 2 - Most DSC objects fall on or below this unity relation. 
These clusters are generally consistent with non-equilibrium 
dynamical conhgurations, and may require at least one other 
effect to explain their high Y 0 yn values, e.g. dark gravitating 
components. On the other hand, a GC with v sin i/<J \/2 v/3 < 1 
can have net angular momentum that can provide a stable 
configuration against rotational breakup. In this case, GC 
rotation alone can account for their Ty yn . Significant error 
bars notwithstanding, all of the “classical” GCs along with 
those on the DGTO sequence are exclusively consistent with 
v sin i/cri/ 2 V / 3 < 1. Thus, their elevated T y' n could be ex¬ 
plained by rotation without the need to invoke additional com¬ 
ponents, bolstering the interpretation that they represent the 
high-luminosity tail of NGC5128’s GCLF. 

As a corollary of the previous exercise, we investigate in 
the right panel of Figure 9 whether the ratio v smi/<Ji 
correlates with .M 1/2 (see Section 3.6). We find that the DSC 
and DGTO GC samples appear to exhibit correlations in this 
parameter space, hinting at rotational support that increases 
with GC mass. This result concurs with the recent findings 
of Kacharov et al. (2014) and Fabricius et al. (2014), who 
measure small but significant rotational speeds in MW GCs. 
However, given that sin i is likely to be randomly distributed 
in the range [ 0 .. 1 ] (hence, (sin i) = 2/7 r), these correlations, 
if real and due to rotational support alone, should be much 
noisier than what we see in the right panel of Figure 9, and 
are therefore probably driven by other effects than rotation 
alone. 


4.3.2. Central Massive Black Holes 


A potential source of artificially enhanced .Mdyn values are 
the effects of central intermediate-mass black holes (IMBHs; 
e.g. Safonova & Shastri 2010; Mieske et al. 2013; Leigh 
et al. 2014) of lesser mass, but otherwise not unlike that 
found recently in a UCD (Seth et al. 2014). To estimate 
the influence of a putative central compact object we com¬ 
pute expected IMBH masses by using our <to estimates with 
the BH mass vs. velocity dispersion relation, Mbh — < 7 , for 
CSSs, which is offset from that of pressure supported galac¬ 
tic systems (Mieske et al. 2013, their Figure 6). From the 
-M[)j[ — o relation, assuming that it scales to lower-mass stel¬ 
lar systems, we obtain BH masses in the following ranges 
6.6 • 10 4 < Mbr/Mq < 5.0 • 10 7 for the DGTO GCs, 
8.5 • 10 2 < Mbh/Mq < 2.2 • 10 6 for the DSC objects, 
and Mbb < 7.6 ■ 10 5 M 0 for the combined “classical” and 
intermediate-Y(V n sample. 

For each GC we integrate the central stellar light pro¬ 
file using our numerical models (see §3.5) until the radius 
of the sphere encompasses a stellar mass corresponding to 
2 -Mbh> assuming the median = 2.2. This radius de¬ 

fines the IMBH sphere of influence (r, ; Merritt 2004). We 
then compute the fraction of the stars within this sphere with 
respect to the modeled population falling within our apertures, 
f(r<Ti). In Figure 10, we plot this fraction as a function of 
-Mdyn for our entire sample. 

While most of the “classical” and intermediate-Yj 1 /" GCs 
show little dynamical influence by potential IMBHs, there is 
an upturn in /(r < r^) for GCs with Mdyn 10 6 M 0 . The 
high-mass sub-sample is mostly made up of the DGTO and 
DSC objects. A Mood test (Mood 1950) for equal medians 
(see horizontal lines in Figure 10) provides no evidence for 
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Fig. 10.— Fraction of the stellar mass within the sphere of influence of a 
putative central IMBH, as a function of A'ldyn for our entire sample. The 
black dashed line shows the median of the combined DGTO and DSC se¬ 
quence GCs, while the red- and blue-dotted lines show the medians for the 
individual sub-samples. 

a difference in the DGTO and DSC sequence f(r < r,) Ip- 
valued.00). Specifically, Figure 10 shows a bimodal distri¬ 
bution in f(r < n) for the DSC sequence objects, as 7/17 
and 8/17 fall either at the ~ 100% or < 20% levels, respec¬ 
tively. Two of these objects show intermediate f(r < ri) ~ 
50 — 60%. For the DGTO sequence GCs there is a smoother 
transition from those with little to no influence by a putative 
IMBH (7/20) to those which would have the majority of their 
stars dynamically dominated by such an object (8/20). Taken 
together, the presence of central IMBHs could in principle 
provide a plausible explanation for the observed dynamics of 
many objects on both sequences. 

We point out that some GCs in the grey shaded region in 
Figure 10 are identified as X-ray sources, almost all of which 
are classified as Low-Mass X-ray Binary (LMXB) hosts from 
Chandra observations (Liu et al. 2011). 2/17 DSC objects are 
X-ray sources, compared to 8/20 DGTO GCs. This is in line 
with the following argument: given the fainter nature of the 
DSC objects, they must have fewer stars compared to DGTO 
clusters at a given ATd yn - With fewer stars providing stel¬ 
lar winds/mass-loss, one would generally expect accretion 
onto an IMBH/LMXB to be less likely compared to DGTO 
sources. 

Given the apparently enigmatic properties of the faint DSC 
subsample, it is important to note that BHs dominating their 
dynamics might alter our basic assumption of the canoni¬ 
cally accepted King (1966) stellar density profile. Spatially 
resolved a profiles would test the putative central BH sphere 
of influence and thus the validity of our assumption. Com¬ 
pared to a uniformly distributed mass component, a central 
IMBH can mimic a dynamical mass as much as 4 — 5 x higher 
than the mass of the BH itself (Mieske et al. 2013). Scaling 
AAdi /2 in this sub-sample down by four then suggests that 
BHs of masses 3.9 ■ 10 4 < ATbh /Mq < 1.3 • 10 6 could plau¬ 
sibly provide the er ppx f that we observe, as well as the elevated 


Y 4 ; yn values. In any case, the lack of a strong correlation in the 
f(r<ri) vs. ALiyii plane calls into question whether IMBHs 
would be the only driver of the two sequences. 

Despite the plausibility of single massive central BHs ex¬ 
plaining some of our results, the potential connection to DSCs 
predicted recently by Baneijee & Kroupa (2011) needs to be 
considered. In this scenario, neutron star and BH remnants of 
massive stars gather as a very concentrated central dark sub¬ 
cluster. Passages through a strong tidal field act to strip lumi¬ 
nous matter from the outskirts, resulting in very high mass- 
to-light ratios. To be observable, the stellar stripping process 
must act on timescales shorter than the self-depletion of dark 
remnants via encounter-driven mechanisms (e.g. three-body 
interactions). Banerjee & Kroupa (2011) predict lifetimes of 
such objects with stellar masses < 10 5 Mq orbiting within 
5 kpc of a MW-like potential to be generally less than 1 Gyr, 
calling into question the likelihood of observing such a cur¬ 
rent population. With that said, we note that given our DSC 
stellar masses of AT* ~ 1 — 5 • 10 5 Mq, and R gc > 5 kpc, 
combined with the predicted correlation between lifetime, and 
both AT* and i? gc , these clusters could plausibly have DSC 
phases lasting on Gyr timescales. Detailed future modeling 
will be critical in determining the plausible parameter space 
necessary for the existence of these objects given the unique 
history of NGC 5128. 

4.3.3. Dark Matter Halos 

If central IMBHs are not driving the dynamics of the DSC 
sequence, then explaining objects with Y 4 /' 1 > 10 becomes 
very difficult without requiring a significant amount of dark 
matter (DM). It is generally accepted that “classical” GCs 
are devoid of non-baryonic matter, but this cannot be entirely 
ruled out (e.g. Moore 1996; Conroy et al. 2011; Sollima et al. 
2012; Ibata et al. 2013) and may in fact be expected from the¬ 
oretical considerations (e.g. Peebles 1984; Saitoh et al. 2006). 
Having considered and ruled out inflated Y 4 /" values due 
to myriad observational/instrumental effects (see Appendix), 
and/or severely out-of-equilibrium dynamical states (see Sec¬ 
tion 4.3.1) due to, e.g., rotation, we now consider the impli¬ 
cations of the DSC sequence being due to the onset of DM 
domination in low-mass systems (AT >10 5 Mq). 

The red-dashed line in Figure 11 shows the difference in 
Yy yn between the DGTO and DSC power-law fits as a func¬ 
tion of AT d yn. In other words, this relation shows the amount 
of DM required to explain the exponentially increasing mass 
excess compared to “classical” GCs. If dark matter is behind 
the DSC objects, then the sharp truncation at ATd yn ~ 10 5 Mq 
hints at a population of increasingly DM dominated structures 
in the immediate vicinity of NGC 5128 with masses as low as 
a few times 10 5 Mq. 

Recent modeling has shown that in a realistic galaxy cluster 
potential, as much as 80-90% of a dwarf galaxy’s DM halo 
may be stripped before baryonic losses become observable 
(Smith et al. 2013a). If the DSC objects originate from low- 
mass halos that have been stripped during their passage(s) 
through NGC5128’s potential well, this then implies pro¬ 
genitors of ATdyn (and thus Y 4 7 n ) ~ 10x higher than we 
observe. This scenario would suggest that at one point in 
time during the assembly history of NGC 5128 there may 
have been a population of compact, low-luminosity baryonic 
structures inside the virial radius of NGC 5128 embedded 
within extraordinarily compact dark matter halos of 10 7 < 
ATdyn/ATg < 10 8 with 100 < Y^/MqLq 1 < 700 (e.g. 
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Fig. 11.— Difference in T77 n between the red- and blue-dashed lines in 
Figure 6 (red relation), and the observational limit of our data (black rela¬ 
tion) as a function of Mdyn- The dominance of the dark gravitating mass 
component exponentially increases above a truncation mass of a few times 
10 s Mq. 

Ricotti 2003). 

If the above is true, then the truncation of the red-dashed 
line in Figure 11 at ~ 10 5 .Mdyn may mark the limit be¬ 
low which primordial DM halos have not survived accretion 
events onto larger galaxy structures. This result is consistent 
with the picture recently put forward in Puzia et al. (2014) 
where the sizes of outer halo star clusters are truncated by 
the abundance of small DM halos (see also Carlberg 2009; 
Carlberg et al. 2011; Carlberg & Grillmair 2013). In this sce¬ 
nario, CSSs of larger masses suffer higher levels of dynamical 
friction than low-mass clusters while moving through the po¬ 
tential well of the host galaxy and sink closer in to the central 
body (e.g. Lotz et al. 2001). While sinking, they suffer extra 
harassing encounters in the more crowded core region, which 
act to further truncate their sizes (e.g. Webb et al. 2013). The 
above may be the explanation for why the DSC objects have 
sizes typical of Milky Way GCs (see Figures 6 and 7, and also 
Figure 16 in Puzia et al. 2014). 

4.4. Potential Progenitors of the DSC Sequence GCs 

Figure 7 shows that with magnitudes of — 8.5 < My < 
—7.4 mag, the DSC objects are among the intrinsically 
faintest star clusters of our sample. This range scratches the 
peak of the GCLF and begins to infringe upon the realm of 
ultra-faint dwarf galaxies (UFDs) in the Local Group, which 
are thought to be the extension of the dwarf spheroidal (dSph) 
galaxy population to lower luminosities {My f —8 mag; e.g., 
Willman et al. 2005; Zucker et al. 2006; Belokurov et al. 2007; 
Zucker et al. 2007; McConnachie et al. 2009; Brown et al. 
2012). Despite smaller sizes than known Local Group dSphs, 
the combination of high Ty 11 and low luminosities shown 
may suggest an evolutionary link to a putative population of 
dSph- and/or UFD-like, dark-matter dominated dwarf galax¬ 
ies. If so, they must have undergone an extremely compact 
and therefore efficient early star-formation burst before being 


tidally stripped of most of their dark-matter halos during sub¬ 
sequent interactions with NGC 5128 (e.g. Smith et al. 2013a; 
Garrison-Kimmel et al. 2014; Sawala et al. 2014). 

The fact that we do not observe current tidal features indica¬ 
tive of stripping might be due to the surface brightness limits 
of available data, but higher (i.e. non-equilibrium) LOSVDs 
might still be consistent with theoretical predictions. For in¬ 
stance, Smith et al. (2013b) model a dwarf galaxy similar 
to the UFD UMall, which with cr = 6.7 km s -1 (implying 
Yy yn ss 1000; Simon & Geha 2007) is an obvious likely can¬ 
didate for a DSC sequence object progenitor. By subjecting 
it to tidal forces in a Milky Way-like gravitational potential, 
these authors successfully reproduced many observed prop¬ 
erties (e.g. luminosity, central surface brightness, ellipticity, 
etc.) of UMall, and found that the galaxy’s cr can be boosted 
on timescales of a few Gyr, in particular around the apocen- 
tre of the orbit. Their modeling showed that “cr boosting” can 
easily reproduce UMall’s cr, and can even reach extreme lev¬ 
els of > 20 km s -1 when the orbital trajectory is close to being 
perpendicular to the disk of the host galaxy. The models there¬ 
fore suggest that the extreme Yy 11 for the DSC objects could 
potentially be explained without requiring large amounts of 
DM. While this result, assuming that “cr boosting” scales to 
a gE galaxy like NGC 5128, might be sufficient to explain a 
handful of the DSC clusters, we consider it highly unlikely 
that all such clusters can be accounted for, due to the required 
synchronization of their apocentre passages and current line- 
of-sight alignments. Note that this mechanism would gener¬ 
ally introduce scatter rather than producing the sequences ob¬ 
served in Figure 6 . 

Assuming that “cr boosting” is unable to solely account for 
the DSC objects, then the small sizes compared to UFD-like 
galaxies call for careful skepticism if DM is the preferred ex¬ 
planation. For example, characteristic central DM densities 
for dwarf galaxies, if canonical halo profiles are assumed, can 
range from 0.1 Mq pc ” 3 for preferred cored profiles, to as 
much as 60 Mq pc ” 3 for cuspy profiles (e.g. Gilmore et al. 
2007; Tollerud et al. 2011). In this case, the required DM 
masses (A4 dm) of 1.6 • 10 5 < Mdm/Mq < 5.4 • 10 6 and 
corresponding densities (pdm) of 40 < pdm/MqPC ” 3 < 
4.5 • 10 4 within the DSC cluster half-light radii are at least 
three orders of magnitude higher than those expected from a 
cuspy dwarf galaxy profile, or if cored, they would each need 
to be embedded in ~ 10 14 — 10 15 Mq halos. 

Irrespective of this problem, such high DM densities might 
actually give rise to central baryonic concentrations that en¬ 
able the formation of IMBHs. This might be realized by trig¬ 
gering an extremely dense central stellar environment, leading 
to a runaway collisional event that culminates with the forma¬ 
tion of a central BH that might continue to grow through the 
accretion of binary star and higher-order multiplets. 

Given all of these interpretations, it seems just as likely 
that a combination of central IMBHs and/or dark stellar rem¬ 
nants, cuspier-than-expected DM halos, and/or “cr boosting” 
could be at work. This conclusion may in fact be the most 
reasonable one, given that IMBHs would presumably be of 
dwarf galaxy origin given the difficulty in building up such a 
mass in a GC-like structure without just as remarkable early 
star-formation efficiency. Then, the residual DM leftover after 
stripping, if cuspy, in combination with “cr boosting”, may be 
sufficient to produce these new objects. This final interpreta¬ 
tion would then require both less extreme BH and DM masses, 
while avoiding the orbital synchronization problem. Regard- 
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less of the mechanism giving rise to them, the emergence of 
the DSC sequence objects is a very unexpected result and calls 
for detailed follow-up observations and high spatial resolution 
modeling. 

5. SUMMARY & FUTURE OUTLOOK 

New high-resolution spectra of compact stellar objects 
around the giant elliptical galaxy NGC5128 (Centaurus A) 
were analyzed. We combined these data with a re-analysis 
of 23 clusters from the literature and used a penalized pixel 
fitting technique to derive new radial velocities for 125 ob¬ 
jects (3 first-time measurements), as well as line-of-sight 
velocity dispersions for 112 targets (89 first-time measure¬ 
ments). Based on these estimates we derived dynamical mass 
and mass-to-light ratio estimates by combining the new kine- 
matical information with structural parameters (mostly ob¬ 
tained from HST imaging) and photometric measurements 
from the literature. 

We briefly summarize our results as follows: 

• At intermediate GC masses (10 5 < A4d yn /Mq < 10 6 ) we 
find the expected population of “classical” GCs, with no 
anomalous kinematical results. These GCs resemble those 
of the Local Group in every way, albeit being slightly 
brighter than average due to our sample selection bias. 

• At the high-mass end (A4d yn ^ 10 6 Mq), we find at least 

two distinct star-cluster populations in the Ty yn -Afdyn 
plane which are well approximated by power-laws of the 
form Ty yn cx Af ^ yn . 

- The “DGTO sequence” is comprised of objects with 
A4 dyn > 2 • 10 6 M & and T^ yn < IOMqL - 1 and 
is well described by a power-law with a slope a = 
0.33 ±0.04. 

- The “DSC sequence” objects have A4d yn similar to 
the DGTO sequence clusters, but with a significantly 
steeper power-law slope, a = 0.79 ± 0.04. More¬ 
over, the faint magnitudes (My > —8.5 mag) lead 
to anomalously high (Yj lyn > 6MqLq X ) mass-to- 
light ratios. We point out that at least two LMC GCs 
(NGC 1754 and NGC2257) also appear to follow this 
relation. 

Despite being among the brightest clusters of our sample, 
some DGTO sequence objects show relatively high mass-to- 
light ratios in the range 5 < T y rn /MqLq 1 < 10, which 
require explanation if a non-baryonic mass component is to 
be avoided. We find that extreme rotation and/or dynamically 
out of equilibrium configurations can explain their kinemat¬ 
ics as well as indications of higher levels of rotational sup¬ 
port with increasing ,V4,i ytl . Plausible alternative and/or con¬ 
current explanations also include particularly top- or bottom- 
heavy IMFs, or the dynamical influence of central IMBHs. 
Altogether we consider that these objects represent the very 
bright tail of the GCLF which is well represented around 
NGC 5128, but poorly populated in the Local Group. 

While the DGTO sequence has a fairly pedestrian expla¬ 
nation, the 6 < T'y’^/MqLq 1 70 values of the DSC se¬ 
quence objects are much more difficult to reconcile with the 
DGTO branch scenarios, save for a small subset. We investi¬ 
gated in detail (see Appendix, for details) the potential impact 


of observational, low S/N, and/or instrumental effects on arti¬ 
ficially inflating the DSC values, and found that astrophysical 
explanations are required. 

For most of these objects the average stellar velocities (vis- 
a-vis the observed op values) do not appear to be high enough 
for extreme rotation to explain their dynamics, and if they 
were, then the clusters would be unstable against rotational 
break-up. Moreover, it is highly unlikely that all clusters 
would have their rotation axes aligned with the plane of the 
sky (assuming sin / £ (7[0,1]), as would be needed to repro¬ 
duce and minimize the scatter in the observed Tp yn vs. A4d yn 
relation. Combined with the difficulty in explaining a mech¬ 
anism that would impart and maintain the necessary angu¬ 
lar momentum for such large rotational velocities, we thus 
consider extreme rotation/significantly out-of-equilibrium dy¬ 
namical configurations insufficient and unlikely to explain 
their properties. 

We considered the plausibility of central IMBHs and a 
central accumulation of dark stellar remnants, consisting of 
stellar-mass BHs and neutron stars, to be driving the extreme 
dynamics of the DSC clusters. If they exist, then putative cen¬ 
tral IMBHs and remnant population can plausibly influence 
enough of the stars in many of the DSC clusters to provide 
an explanation for their velocity dispersion measurements. In 
fact, if our assumed canonical structural profiles were suffi¬ 
ciently perturbed by an IMBH’s presence and/or stellar rem¬ 
nant population, then this explanation requires even less mas¬ 
sive IMBHs to become plausible. With that said, given that 
the IMBH + stellar remnant interpretation can only account 
for the dynamics of ~ 50% of these objects, this is unlikely 
to explain the emergence of the DSC sequence. 

If central IMBHs and stellar remnant populations are not 
the only cause of the elevated Yy yn , then the possibility of 
significant dark matter mass components must be considered, 
despite the wide acceptance that “classical” GCs are devoid 
of DM. This result would have important implications for GC 
formation models and early structure formation, and indicate 
that not all extragalactic star clusters are genuine GCs. More 
importantly, the presence of such amounts of DM in the DSC 
sequence objects would imply that they represent the lowest 
mass primordial dark matter halos that have survived accre¬ 
tion onto larger-scale structures to the present day. In other 
words, there may still exist a large reservoir of 1CP Mq- scale 
dark matter halos surviving in relative isolation in the universe 
today, at least around relatively quiescent larger dark matter 
halos like NGC 5128. Moreover, if these objects are stripped 
of formerly more massive dark matter halos, presumably as 
former dwarf galaxies, this would imply the presence of a sig¬ 
nificant collection of objects with 10 7 < A4d yn /A/o (S 10®, 

and 100 < T ( y n /M qL~^ 1 < 1000 in the relatively recent 
past of NGC 5128. 

This interpretation is not without its serious problems. For 
example, with the above properties, central (~ 10 pc) DM 
masses/densities on the order of ~ 10 3 x larger than canonical 
DM halo profiles would be required; a scenario that cannot 
be reconciled with any current theoretical framework. Given 
the improbability that such massive central BHs, exotic and 
ultra-concentrated DM halos, or extremely out-of-equilibrium 
dynamical configurations can individually explain the proper¬ 
ties of the DSC sequence objects, it seems perhaps most likely 
that a mixed bag of such factors may be at play behind this re¬ 
sult, although it is puzzling how a combination of these phys¬ 
ical mechanisms would conspire to generate a relatively sharp 
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Yy yn vs. A^dyn relation. 

It remains to be seen if similar objects exist in the star clus¬ 
ter systems of other galaxies, but verification of these results 
will be difficult for more distant systems due to their intrinsic 
low luminosities. Nonetheless, detailed chemical abundance 
studies of these objects will shed light on the origins (e.g. pri¬ 
mordial or not, simple or multi-generational stellar popula¬ 
tions, etc.) of both DGTO and DSC clusters. While the prox¬ 
imity of NGC5128 provides the possibility to study the in¬ 
ternal dynamics and stellar populations of its CSSs, the dis¬ 
tance approaches the limits of what is currently feasible with 
today’s instrumentation on 8-10m class telescopes. Nonethe¬ 
less, large-scale, complete studies of the chemo-dynamics 
of GCs in the Local Group, and of giant galaxies within 
~ 5 Mpc, will help reveal the true nature of this enigmatic 
new type of compact stellar systems. 
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APPENDIX 

TESTING FOR POTENTIAL DATA ANALYSIS BIASES 

We investigated the potential of erroneous literature measurements giving rise to the DSC sequence in Figure 6 , and find that 
one object may be explained by discrepancies found in the literature. For GC0225, we took the apparent magnitude V = 19.93 
mag from the Woodley et al. (2007) catalogue, which is 2.86 magnitudes fainter than in the discovery publication of Holland et 
al. (1999). We do not attempt to explain this discrepancy, but note that if the brighter measurement is used, it leads to Y p '" a 
factor of ~ 14 lower, bringing it more in line with the DGTO sequence. Additionally, Mieske et al. (2013) note an inconsistency 
in the size measured by Holland et al. (1999), suggesting that 77 , is ~ 3.8x smaller than originally estimated. The smaller size 
would naturally lower Ty yn by the same factor, bringing it more in line with the two intermediate (but still > 10) Y p yn clusters. If 
GC 0225 is both brighter and smaller, then its Tj 1 ; n would approach one. Eliminating it from the power-law fit results in virtually 
the same slope (ct = 0.78±0.06), with slightly larger scatter. 

We also tested for various data reduction effects that may artificially give rise to the DGTO and DSC sequences. For example, 
a straight-forward explanation is that in performing the convolutions to estimate <T ppx f, the ppxf code may have “jumped" over 
the targeted Mg 6 and Fe 5270 spectral features and based the kinematics on the wrong combination of spectral lines due to the 
relatively low-S/N of some spectra. This explanation would require an error on u r ppx f of > 100 km s -1 which Figure 4 indicates 
is not present in our data. Still, given that multiple v r measurements can be found in the literature for many of our sample GCs, 
we individually investigated the DSC sequence objects to search for any literature v r estimates that are discrepant by > 100 km 
s -1 . We found only three such clusters (GC 0115, GC 0225, and GC 0437 by Ai> r , max =142, 93, and 110 km s” 1 , respectively), 
noting that they each fall within the error bars of at least one literature value. Nonetheless, even if all of three GCs are omitted 
from the fit shown in Figure 6 , the relation again does not change significantly (a = 0.78±0.06). 

TESTING FOR TARGET CONFUSION 

The significant luminosity difference between the DGTO and DSC sequence GCs calls for investigations into whether the 
bifurcation is due to observational effects. In general, it is unlikely that any contamination of the sampled GC flux by other 
sources would produce two such relations, instead of just increased scatter. However, to address this issue we have visually 
inspected all available archival HST imaging data for any potential contamination of the DSC sequence GCs and found no 
indications for any target confusion due to foreground starlight, background galaxies, or enhanced surface brightness fluctuations 
of the surrounding galaxy light for each GC. 

TESTING FOR SPURIOUS RESULTS DUE TO NOISY SPECTRA 

The well-defined DSC relation could be due to the observational limit of our data. The solid black line shown in Figure 6 shows 
where GCs with My = —7.4 mag (the faintest measured objects in our sample) would lie. Objects to the left are inaccessible to 
our survey, and the nearly parallel alignment of this line to the DSC sequence suggests that the small scatter may be an artifact of 
this limit. Indeed, the left panel of Figure 12 shows cr ppx f as a function of My for our entire sample with symbols as in Figure 6 , 
but GCs with unreliable cr ppx f shown as open grey circles, and the inset histogram showing the My distribution of GCs for 
which ppxf could not derive a cr ppx f estimate. It can be seen that most objects with poor/unavailable <r ppx f encroach upon the 
luminosities of the DSC sequence. These objects are possibly “classical" GCs with absorption features too narrow for GIRAFFE 
to resolve, and/or for ppxf to accurately measure through the noise, whereas DSC objects have features sufficiently broadened to 
be measurable. This effect would give the false impression that almost all objects of My > —8.0 mag seem to have anomalously 
high Yj 1 '". In fact, the objects on the DSC sequence (red squares) show a trend of higher luminosity with larger <r ppx f, contrary 
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Fig. 12.— (Left panel): cr ppx f values as a function of My , with brightness increasing to the right. Symbol shapes are as in Figure 6, with the addition of 
objects for which <7 ppx f could not be estimated (grey inset histogram) or which we considered unreliable (open grey circles). (Right panel): Difference between 
velocity dispersion measured on spectra with artificially added noise, cr no i sy , and our measurements, cr ppx f , as a function of the amount of noise added in units 
of the standard deviation of the original spectra. The errors of the cr noisy measurements are shown on a logarithmic scale, parametrized by the color shading 
to illustrate the increased uncertainties as more noise is added. If poor-quality spectra were the cause of the DSC sequence, an upturn of the c7 no i S y — cr ppx f 
difference would be expected towards larger noise factors; however, the symmetric increase in spread to both negative and positive cr no i sy — cr ppx f values clearly 
rules out this possibility. 


to the expectation if noise were “washing-out” the finer spectral details used to estimate <T ppx f. 

Since the DSC sequence objects are among the faintest objects in our sample (see Figure 7), an obvious point of concern 
is that their spectra are among the noisiest. We therefore performed the following test to check against the potential of noise 
“washing out” some of the finer spectral details used to estimate er ppx f, thus tricking the code into estimating systematically 
wider dispersions. We took four DGTO sequence GCs (GC 0106, GC 0277, GC 0306 and GC 0310; chosen to cover a reasonable 
luminosity range), artificially added varying levels of noise, and repeated the er ppx f measurements. For each GC we added noise, 
pixel by pixel, by randomly drawing from a N(0,n • cr spec ) distribution where cr spec is the flux dispersion intrinsic to each 
spectrum, and n is drawn from 17(0, 5) distributions. This procedure decreases the spectral S/N by a factor 1/n. Repeating this 
process 100 times for each GC then builds a picture of how well behaved the ppxf code is for increasingly noisy spectra. The right 
panel of Figure 12 shows the results of this exercise for all four GCs and it is clear that even for large amounts of noise, the points 
cluster well around the adopted <r ppx f values, with increased dispersion bracketing the nominal values accompanied by larger 

errors. If sudden “jumps” in er ppx f due to noisy spectra were the cause of the elevated Ty 11 shown by the DSC sequence, then it 
would be expected that the points in Figure 12 would tend towards higher (T no i sy estimates with stronger amplified noise. No hint 
of such asymmetric trend is seen in the plot, and we conclude that the DSC sequence is not artificially created by poor quality 
spectra. 


TESTING FOR CORRELATIONS WITH GALACTOCENTRIC RADIUS AND AZIMUTHAL ANGLE 

In searching for the origins of the structures found in Figure 6, in particular the bifurcation at high Aidyn, we test whether any 
correlation of the measured er ppx f with spatial distribution might give deeper insight into our results. For the following, we note 
that we are using the raw <r ppx f measurement directly obtained from the spectra without subsequent correction for GC light profile 
sampling. We plot in Figure 13 er ppx f as functions of galactocentric radius ( R gc ; left panel), and azimuthal angle (<f>; right panel), 
where 0° corresponds to north and $ increases clockwise. We find no indications for any er ppx f--R gc or cr ppx f-'f > correlations for 
our GC sample. More importantly, we find no such correlations for either of the GC sub-samples, in particular for the DSC GCs 
(red squares). We conclude that systematic variations in <r ppx f as a function of spatial distribution are not responsible for the 

bifurcation in the Ty yn -Addyn relation. 

TESTING FOR GALAXY SURFACE BRIGHTNESS CONTAMINATION 

Here we investigate the fractional flux from NGC5128’s diffuse galaxy light component that enters each fibre together with 
the GC light. We consider the potential impact of such fibre contamination, which could potentially lead the ppxf code to predict 
artificially high er ppx f values. However, if this was the case, one would expect that the associated error bars would reflect this 
uncertainty (given full error propagation), and/or that the elevated Ty yn values required to explain the DSC sequence (see §4.2) 
would either tend towards higher values with decreasing luminosity, which is not seen in our data (see Figure 7 and 11), or that 
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Fig. 13.— (Left panel)'. Raw cr ppx f measurement shown as a function of galactocentric radius, Rgc, with symbols as in Figures 6 and 7. (Right panel): Same 
as in the left panel, but showing cr ppx f as a function of azimuthal angle, with 0° corresponding to North and increasing clockwise (see also Figure 1). No obvious 
correlations between cr ppx f, Rgc, or $ are seen, effectively ruling out the possibility of systematic variations of <r ppx f with spatial distribution giving rise to the 
DGTO or DSC sequences. 



the right panel with the decadal logarithm of the fractional NGC 5128 galaxy light contribution to each of the sampled GC fluxes, before sky subtraction. 

there would be a vanishing or no T' t l ''"-Ad,i yn correlation at all. As can be seen in Figure 6, the DSC sequence is consistent with 
the individual error bars, which themselves are dominated by the fully propagated uncertainties on a. 

Given that our background light SED modeling and subsequent subtraction from each fibre is not a localized process (see 
§ 2.3), visual inspection of the sky+background subtracted GC spectra indicated that they should not be significantly affected 
by residual galaxy light contamination. To be sure, we explicitly calculated the expected residuals by comparing the difference 
between the flux of NGC 5128 at the location of each target and our final sky estimates to the flux of our final target spectra. The 
results suggested that residual components due to diffuse galaxy light is for all GCs < 4%, and in most cases near the 1% level, 
with the expected trend toward lower residuals with increasing R gc . We then tested for the maximal effects of these residuals 
on the DSC sequence objects by adding the signal from the nearest sky fibre to the final, reduced target spectra at flux levels 
between 1% and 90% of the object flux and re-measuring er ppx f. The results showed that below ~ 20 %, the effect on cr ppx f was 
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negligible, above which <r ppx f monotonically growing with increased sky contamination, as one would expect. 

In any case, we plot in the left panel of Figure 14 the Ty yn -A / ldyn plane, as in Figure 6, but parametrize the data point 
color this time with the local surface brightness of NGC5128 at each GC location. We obtain the surface brightness values 
from the profiles measured as part of the Carnegie-Irvine Galaxy Survey (CGS; see Ho et al. 2011; Li et al. 2011). No obvious 
correlations with local surface brightness values are found that would drive the bifurcation in the Y p y " -A4 d yn relation. We 
continue by using the surface brightness profile for NGC5128 to compute the fractional background flux, /(/,(>„a entering 
each fibre with the GC fluxes. The right panel of Figure 14 shows the results, with color parameterized by log /(LcenA) before 
any sky+background subtraction is performed. Several results are shown: 1) The “classical” and DGTO sequence GCs reveal 
the expected correlation between GC mass (i.e. luminosity) and /(LcenA) with a decreasing /(LcenA) sequence from ~ 30% 
for the low-mass GCs to negligible fractions for the highest-mass clusters; 2) GCs on the DSC sequence show on average larger 
f(LcenA) than GCs on the DGTO sequence as expected from the trends seen in Figure 7, which shows that DSC GCs have 
luminosities of —8.5 < My < —7.5 mag; 3) The DSC sequence GCs show no correlation between / (Lcen a ) and Yy yn or 
Miy-a, in contrast to the expectation if the T p , yn -A / ldyn sequences were created by insufficient background light subtraction. 

While Figure 14 illustrates fibre contamination before any sky modeling and subtraction was carried out, we performed 
one final test to check for the worst effect that insufficient background subtraction could have. We first identified higher- 
S/N GCs with CT ppx f equivalent to that which would give a “normal” T P “ for a DSC sequence cluster. For example, given 
the DSC GC0324’s cr ppx f,GC0324 = 24.89 kms -1 and Ty y g C0324 = 14.07 MqLq 1 , we identified the DGTO GC0050 with 

Cppxf,GC 0050 — 10.09kms -1 , which would give Y' \sqco 324 = 2.31 MqLq 1 . Using three such combinations, we scaled the 
individual spectra of the bright GC, pre-sky subtraction, to the flux level of the fainter cluster, and added the nearest sky fibre’s 
signal. Performing the sky subtraction as before (see § 2.3) with the next three closest fibres then simulated the noise levels of the 
faint GC, but preserved the spectral details of the brighter. Measuring cr ppx f on the degraded spectrum was, at worst, consistent 
with our adopted measurements, and even led to a slightly lower cr ppx f for two of the experiments. We conclude that improper 
background subtraction is thus very unlikely to inflate our cr ppx f measurements of the DSC sequence objects. 
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